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Preface 


If  a  television  commercial  by  General  Electric,  which  aired  in  May  2003,  about  a 
nanotechnologist  and  supermodel  who  fall  in  love  at  first  sight  were  to  be  taken  as 
evidence,  we  could  conclude  that  nanotechnology  has  replaced  rocket  science  as  the 
discipline  of  the  ultra-chic  uber-geek.  Movies  such  as  Spiderman  and  Minority  Report, 
and  mass  market  literature  like  Prey  have  shaped  public  consciousness  of  the  term  so  that 
it  elicits  images  of  tiny  miracle  working  robots  and  lethal  gray  goo.  In  fact,  Merriam- 
Webster’s  Online  Dictionary  gives  this  definition:  “the  art  of  manipulating  materials  on 
an  atomic  or  molecular  scale  especially  to  build  microscopic  devices  (as  robots).”  But 
most  serious  research  in  the  field  of  nanotechnology  is  not  about  trying  to  build  tiny 
robots,  and  there  are  plenty  of  less  glamorous  disciplines  like  chemistry  or  materials 
science  that  are  focused  on  manipulating  materials  on  an  atomic  or  molecular  scale.  So, 
what  separates  nanotechnology  from  other  buzzwords  like  “lockbox”  or  “smoking  gun” 
which  are  just  plagiarisms  of  emotionally  charged  language  with  little  substance? 

For  me,  nanotechnology  is  meaningful  and  different  from  other  scientific 
disciplines  because  it  is  about  engineering  solutions  to  problems  that  take  advantage  of 
the  unique,  controllable  properties  of  nanoscale  matter.  When  the  characteristic  length 
scale  of  an  object  is  a  billionth  of  a  meter,  quantum  mechanics  allows  for  physical 
properties  that  do  not  exist  in  the  macroscopic  world.  For  the  first  time  in  history,  we  can 
change  some  parameter  of  a  simple  nanoscale  system,  which  alters  some  system  property 
in  a  well  understood  way,  so  that  we  can  achieve  some  a  priori  design  goal. 
Nanotechnology  is  a  notoriously  difficult  word  to  define,  and  I  am  certain  that  my 
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definition  is  little  better  than  Merriam-Webster’s,  but  this  is  how  I  differentiate  the  field 
from  others. 

This  thesis  is  about  how  a  particular  nanoscale  particle,  the  nanoshell,  can  be  used 
in  a  system  to  achieve  optimal  energy  conversion  from  electromagnetic  fonn  to  thermal 
or  electronic  form.  I  wrote  computer  code  that  calculates  the  optimal  mixture  of 
nanoshell  geometries,  or  species,  for  absorbing  or  scattering  a  given  spectrum  of  light. 
The  code  is  intended  to  be  sufficiently  general  so  that  it  will  be  a  useful  design  tool  for  a 
wide  variety  of  engineering  problems,  but  I  solved  only  two  very  specific  scenarios 
which  are  of  interest  to  our  work  in  the  Laboratory  for  Nanophotonics  at  Rice  University. 
I  hope  this  work  will  be  valuable  in  the  design  and  construction  of  novel  devices  with 
practical  use,  and  I  hope  it  will  stand  as  a  realistic  example  of  the  current  state  of 
nanotechnology  and  its  near-tenn  potential. 
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Chapter  1  - 

Plasmonic  Particles,  Sunlight,  and  Optimization 


The  Sun  is  the  energy  source  of  life,  and  nature  stockpiled  that  energy  into 
convenient  forms  like  coal  and  oil  over  the  course  of  millennia.  We  are  entering  an  age 
when  depletion  of  these  reserves  causes  us  to  question  the  long  term  sustainability  of  life 
as  we  know  it.  It  is  critical  that  we  develop  sustainable  processes  for  the  production  of 
food,  potable  water,  and  the  extraction  of  resources  required  for  manufacturing  and 
construction.  Success  in  these  goals  hinges  on  our  ability  to  generate  renewable  energy. 
Energy  from  the  Sun  is  abundant,  clean,  and  most  importantly,  sustainable,  but  we  lack 
methods  to  harness  it  efficiently.  New  discoveries  in  nanotechnology  hold  great  promise 
toward  that  end  ’  ’  .  In  this  work,  I  examine  theoretically  how  one  nanoscale  technology, 
the  nanoshell,  can  be  used  to  optimally  absorb  or  scatter  the  energy  in  sunlight,  which 
would  allow  us  to  convert  it  into  usable  electronic  or  thermal  form. 

Nanoshell  Physics 

The  photophysics  of  nanoshells  is  intimately  related  to  the  plasmon  resonance  of 
the  particles.  Plasmons  are  collective  oscillations  of  the  conduction  electrons  in  a  metal, 
and  they  can  dominate  the  metal’s  optical  properties  when  certain  conditions  are  met. 

This  strong  interaction  between  light  and  matter  enables  efficient  coupling  of 
electromagnetic  energy  into  the  metal,  and  the  strength  of  the  interaction  is  determined  by 
the  geometry  of  the  object  as  well  as  its  dielectric  function.  The  dielectric  function 
describes  how  an  electric  field  causes  a  polarization  response4  in  a  given  material.  It  is  a 
measure  of  how  charged  particles  like  electrons  respond  to  light,  so  it  depends  on  the 
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density  of  free  electrons  present  and  the  electron  mean  free  path.  Generally,  the  dielectric 
function  is  complex-valued,  and  it  varies  with  the  frequency  of  the  applied  electric  field. 

The  origin  of  plasmons  is  most  easily  understood  in  terms  of  Drude’s  simple 
classical  description  of  conductivity  in  metal5.  Drude  applied  a  kinetic  theory  usually 
associated  with  the  description  of  an  ideal  gas  to  explain  macroscopic  phenomena  in 
metals.  He  assumed  that  electrons  could  be  treated  as  a  gas  that  is  essentially  free,  but 
that  is  confined  by  the  boundaries  of  the  metal.  The  electrons  are  independent  in  that 
they  do  not  interact  with  one  another,  but  they  do  scatter  off  of  the  positive  ion  cores  of 
the  crystal  lattice.  Drude’s  model  explains  AC  and  DC  conductivity  in  metals,  the  Hall 
effect,  and  thermal  conductivity  surprisingly  well  considering  the  crude  approximations 
he  used. 

The  model  predicts  a  zero  in  the  dielectric  function  at  a  certain  frequency, 

0)B  = 

known  as  the  bulk  plasma  frequency,  where  n  is  the  electron  density,  e  is  the  electronic 
charge,  and  m  is  the  electron  mass,  all  in  CGS  units.  This  frequency  marks  the  boundary 
between  solutions  to  the  wave  equation  that  can  propagate  ( co  >  coB  ,  the  material  is 
transparent),  and  those  that  cannot  ( co  <  coB  ,  the  field  decays  exponentially)  5. 
Furthermore,  when  co  -  coB ,  a  solution  to  a  different  wave  equation,  for  charge  density 
waves,  exists.  At  that  exact  frequency,  collective  oscillations  of  free  electrons,  or 
plasmons,  are  possible. 

More  precisely,  co  -  coB  is  the  condition  that  must  be  met  for  the  electron  charge 
density  to  have  oscillatory  time  dependence.  The  oscillations  arise  because  collective 
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displacement  of  electrons  in  one  direction  will  uncover  the  fixed  positive  ion  cores  of  the 
crystal  lattice.  These  cores  provide  a  restoring  force  on  the  displaced  electron  “gas” 
through  Coulombic  attraction,  and  a  harmonic  oscillation  results.  The  alkali  metals  are 
particularly  well  described  by  the  Drude  model5.  The  noble  metals  are  as  well,  but  to  a 
lesser  extent  because  interband  transitions  from  the  d-  to  sp-orbital  in  the  atomic  lattice 
dominate  the  electronic  response  at  short  visible  wavelengths6.  Other  materials  require 
quantum  mechanical  models  to  explain  their  observed  dielectric  functions. 

Since  plasmons  originate  from  the  Drude  response  of  a  material,  they  are  only 
observable  in  metals  that  have  dielectric  functions  dominated  by  the  Drude  term. 
However,  bulk  plasmons  have  never  been  directly  observed  in  practice,  even  for  alkali  or 
noble  metals,  although  they  undoubtedly  exist.  Bulk  plasmons  are  so  difficult  to  observe 
because  they  cannot  be  directly  excited  by  light.  At  optical  frequencies,  bulk  metals  have 
a  skin  depth  on  the  order  of  a  few  tens  of  nanometers6.  Electrons  on  the  interior  of  the 
bulk  are  left  undisturbed  because  of  the  screening  effect  of  the  surface  electrons.  The 
most  direct  evidence  available  for  the  existence  of  bulk  plasmons  is  the  observation  of 
energy  losses  in  multiples  of  ficoB  when  electrons  are  fired  through  thin  metallic  films7. 

As  mentioned  above,  geometry  also  plays  a  significant  role  in  the  interaction  of 
light  with  matter.  The  easiest  way  to  understand  the  influence  of  shape  is  through 
Maxwell’s  electromagnetic  theory.  Maxwell’s  equations  tell  us  that  the  tangential 
component  of  a  magnetic  field  and  the  normal  component  of  the  electric  flux  density  will 

o 

be  discontinuous  at  a  boundary  between  two  dissimilar  dielectric  media  .  These 
boundary  conditions  give  rise  to  the  importance  of  geometry  in  the  solutions  of  wave 
equations  and  therefore  the  optical  properties  of  materials.  In  the  simplest  case, 
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Maxwell’s  equations  predict  the  existence  of  charge  density  waves  with  a  frequency  of 


along  a  surface  at  the  boundary  of  an  infinite  conductor  in  contact  with  a  dielectric 

medium.  These  waves  are  known  as  surface  plasmons. 

Unlike  bulk  plasmons,  surface  plasmons  can  be  excited  by  light,  but  alas,  not 
directly.  They  require  the  presence  of  an  optical  coupler  that  can  add  enough  momentum 
to  what  is  already  present  in  light  to  excite  traveling  wave  plasmons.  Examples  of  optical 
couplers  include  prisms  or  gratings.  More  information  about  how  these  add  to  the 
momentum  of  light  so  that  surface  plasmons  can  be  excited  is  available  in  Surface 
Plasmons  on  Smooth  and  Rough  Surfaces  by  Raether9. 

Fortunately,  when  a  particle  is  very  small,  comparable  to  the  wavelength  of  light, 
standing  wave  solutions  to  the  surface  charge  density  wave  equation  appear.  These  local 
charge  density  fluctuations  can  be  directly  excited  by  light  and  are  the  phenomena  I  refer 
to  when  I  use  the  word  “plasmon”  hereinafter.  The  fact  that  plasma  oscillations  in 
macroscopic  media  are  so  difficult  to  excite  is  primarily  responsible  for  the  late  arrival  of 
technology  to  exploit  them.  The  advent  of  nanoscale  control  of  metallic  particle 
geometries  enables  scientists  and  engineers  to  use  plasmons  in  a  wide  range  of 
applications  from  medicine10,  to  sensor  technologies11,  to  light  harvesting  and  energy 
conversion12. 

Nanoshells  are  particularly  interesting  plasmonic  particles  because  their  plasmon 
resonance  is  controllably  tunable  over  the  entire  visible  light  spectrum  as  well  as  the  near 
to  mid-infrared  domain.  Nanoshells  are  spherical  core-shell  particles  typically  made 
using  silica  cores  coated  with  thin  gold  or  silver  shells.  In  Figure  1.1,  ri  is  the  radius  of 
the  inner  core,  and  r2  is  the  overall  radius  of  the  particle.  Their  optical  properties  stem 
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from  the  interaction,  or  hybridization,  between  the  plasmon  resonances  of  two  elementary 
geometries,  that  of  a  nanoscale  conducting  sphere  in  an  infinite  dielectric  medium  and  of 
a  nanoscale  dielectric  cavity  in  an  infinite  conducting  medium. 


Figure  1 . 1  The  geometry  of  a  nanoshell 
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Mie  solved  Maxwell’s  equations  for  spherical  geometries  in  1908  .  The  theory 
was  extended  by  Aden  and  Kerker  for  radially  symmetric  core-shell  geometries  in 


1951 l4.  Mie  theory  predicts  that  a  conducting  sphere  will  have  plasmon  resonances  at 


(oB 


2/  + 1 


,  where  /  is  a  positive  integer  quantum  number.  Similarly,  it  predicts  plasmon 


resonances  at  coR 


/  +  1 
2/  +  1 


for  a  cavity.  The  /  =  1  mode  is  known  as  the  dipole  resonance, 


/  =  2  is  the  quadrupole,  and  so  on.  In  a  nanoshell  geometry,  the  elementary  sphere  and 
cavity  modes  hybridize  in  rigorous  analogy  with  molecular  orbital  theory  15.  Plasmon 
hybridization  theory  is  the  best  way  to  gain  an  intuitive  understanding  of  how  to  engineer 
the  optical  behavior  of  nanoshells.  Coulombic  forces  mediate  the  interaction  between  the 
modes,  so  that  the  thickness  of  the  metallic  shell  layer  controls  the  strength  of  the 
hybridization. 

Hybridization  between  sphere  and  cavity  dipole  modes  is  shown  in  Figure  1.2. 


The  cavity  mode  has  energy  equal  to  coB,  1.41  times  higher  than  the  energy  of  the 
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sphere  plasmon  at  -yj^coB  .  In  the  thick  shell  limit,  the  “bonding”  and  “antibonding” 

hybridized  modes  are  only  weakly  coupled,  so  they  are  nearly  equal  in  energy  to  their 
constituent  parts.  The  antibonding  mode  is  known  as  the  dark  plasmon  because  it  has  no 
net  dipole  moment  and  cannot  normally  be  excited  by  light,  while  the  bonding  mode  is 
called  the  bright  plasmon.  As  the  shell  becomes  thinner,  the  Coulombic  force  increases, 
and  the  dark  plasmon  blue-shifts  in  energy,  while  the  bright  plasmon  red-shifts. 

nanoshell  «+ 


Figure  1.2  Hybridization  of  sphere  and  cavity 
plasmon  modes  to  make  the  nanoshell 
plasmon 


In  addition  to  plasmon  hybridization  theory,  resonance  energy  is  affected  by 
overall  particle  size.  Just  as  the  pitch  of  a  large  ringing  bell  is  lower  (longer  wavelength) 
than  smaller  bells,  the  standing  wave  resonance  of  large  nanoparticles  occurs  at  longer 
wavelengths.  The  combination  of  this  effect  with  plasmon  hybridization  enables  the 
unique  tunability  of  nanoshells. 

Since  plasmons  are  excitable  by  light  directly,  their  existence  can  be 
experimentally  verified  by  measuring  light  extinction  in  a  spectrophotometer.  Extinction 
is  the  amount  of  light  scattered  or  absorbed  by  the  particle.  At  wavelengths  near  the 
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resonance  energy,  light  interacts  very  strongly  with  the  particle,  and  the  extinction 
spectrum  will  show  a  pronounced  peak.  Spectra  measured  in  this  way  are  well-fit  by 
extinction  efficiency  spectra  calculated  using  Mie  theory16.  Efficiency  is  the  ratio  of  the 
particle’s  effective  cross  section  to  its  geometrical  cross  section.  Figure  1.3  shows  how 
spectra  calculated  using  Mie  theory  can  be  explained  in  terms  of  the  plasmon 
hybridization  model. 


Wavelength  (nm) 


Figure  1.3  Hybridization  controls  plasmon  resonance  wavelength 
On  the  left,  thicker  shells  result  in  less  interaction 
On  the  right,  thinner  shells  result  in  more  interaction 


The  extinction  efficiency  of  plasmonic  particles  can  range  from  several  times 
greater  than  one  when  on-resonance  to  less  than  one.  From  a  geometrical  optics 
standpoint,  it  does  not  make  sense  for  any  object  to  have  an  effective  cross  section  larger 
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than  its  area  projected  onto  the  plane  of  incidence.  However,  geometrical  optics  breaks 
down  for  nanoscale  objects.  Ultimately,  light  is  just  a  wave  of  time-varying  electric  and 
magnetic  fields,  and  an  object  in  the  light’s  path  interacts  with  all  the  nearby  fields,  not 
just  those  directly  incident  upon  it.  Consider  Figure  1.4,  the  Poynting  vector  diagram  for 
a  small  metallic  particle  near  its  plasmon  resonance.  The  Poynting  vector  is  found  by 
taking  the  cross  product  of  the  electric  and  magnetic  field  vectors,  and  it  represents  the 
flow  of  energy  in  the  system.  The  particle  seems  to  pull  the  light  toward  it  from  all 
directions17, 18.  For  off-resonance  particles,  the  opposite  is  observed.  The  effective  cross 
section  is  really  a  simplification  of  the  complicated  way  light  interacts  with  nearby 
matter. 


Figure  1.4  Poynting  vector  diagram  of  a  plasmonic  particle  on- 
resonance  (left)  and  off-resonance  (right)  (Reprinted  with 
permission.  Copyright  1983,  American  Association  of  Physics 
Teachers17) 


So,  plasmon  hybridization  and  particle  size  explain  the  position  of  the  resonance 
peaks  shown  in  a  plot  of  extinction  efficiency,  and  consideration  of  energy  flow  around 
the  particle  can  explain  the  peak  height.  Only  the  breadth  of  the  peak  remains  to  be 
explained.  Two  factors  influence  peak  width.  The  first  is  the  plasmon  dephasing  time,  or 
the  length  of  time  a  collective  oscillation  of  electrons  can  exist  before  it  disappears  as  a 
result  of  frictional  damping  against  scattering  centers.  Plasmon  dephasing  is  a 
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manifestation  of  the  fact  that  perfect  single-frequency  signals  are  impossible  in  real  life; 
such  a  signal  would  have  to  last  an  infinitely  long  time.  Since  the  plasmon  has  a  finite 
lifetime,  taking  the  Fourier  transfonn  will  result  in  a  peak  that  is  not  perfectly  sharp, 
where  the  peak  width  is  an  indication  of  the  plasmon  dephasing  time.  Second,  in 
ensemble  measurements  plasmon  line  width  is  a  result  of  polydispersity  in  the  nanoshell 
suspension.  There  will  always  be  some  distribution  of  core  size  and  shell  thickness. 
Typically,  cores  might  be  ±5%  from  the  mean  size  and  shells  might  be  ±1%.  An 
ensemble  measurement  of  the  light  extinction  for  such  a  group  of  particles  will  show  an 
average  of  the  many  slightly  varied  plasmon  resonances. 

The  Solar  Spectrum 

In  space,  our  Sun  has  the  emission  spectrum  one  would  expect  from  a 
5800  K  blackbody19'20.  However,  the  light  that  reaches  the  surface  of  the  Earth  is 
scattered  and  refracted  by  the  atmosphere,  so  its  path  from  space  to  the  ground  plays  a 
crucial  role  in  the  spectrum  an  observer  sees.  Since  the  spectrum  on  the  ground  can  vary 
greatly,  the  American  Society  for  Testing  and  Materials  (ASTM)  developed  the  ASTM 
G- 173-03  standard  solar  spectra.  The  spectral  data  can  be  downloaded  for  free  from  the 
National  Renewable  Energy  Laboratory’s  (NREL)  website-  .  Figure  1.5  shows  the 
spectra  with  the  blackbody  curve  overlaid.  Manufacturers  of  photovoltaic  devices  and 
other  solar  energy  technologies  report  perfonnance  in  tenns  of  these  spectra  in  order  to 
facilitate  meaningful  comparisons  between  technologies  or  products.  Usually,  the 
AM  1 .5  spectrum  is  used  because  it  represents  time-averaged  irradiation  over  one  day. 
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2.5 


Wavelength  [nm| 

Figure  1.5  Standard  solar  spectra,  with 
5800  K  blackbody  overlay 

Three  spectra  are  included  in  the  standard:  AM  0,  total  global  AM  1.5,  and  direct 
normal  AM  1.5.  “AM”  stands  for  air  mass,  and  it  is  a  measure  of  the  amount  of 
atmosphere  light  must  travel  through  before  reaching  the  surface.  AM  0  refers  to  the 
spectrum  incident  on  Earth  in  outer  space  before  passing  through  any  atmosphere.  AM  1 
is  the  spectrum  on  the  ground  at  high  noon,  when  the  sun  is  at  its  zenith.  As  the  time  of 
day  continues  past  noon,  the  Sun  will  rotate  away  from  zenith,  and  the  atmospheric  path 
length  of  its  light  increases.  When  the  path  length  is  equal  to  1 .5  times  the  length  at 
zenith,  an  observer  would  see  the  AM  1.5  spectrum.  This  occurs  when  the  sun  is  48.2° 
past  zenith,  as  can  be  seen  from  Figure  1.6.  In  the  figure,  the  curved  atmosphere  is 
approximated  by  a  line  tangent  to  the  outer  atmosphere  and  normal  to  the  AM  1  path, 
forming  a  right  triangle.  As  long  as  the  angle  from  zenith  is  <60°,  the  approximation  is 
valid.  The  difference  between  the  total  global  and  the  direct  nonnal  spectra  is  that  the 
direct  normal  spectrum  only  includes  light  that  arrives  directly  from  the  Sun.  Light  that 
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would  be  incident  elsewhere,  but  is  scattered  toward  the  observer  is  not  included.  The 
total  global  spectrum  includes  all  contributions,  and  it  is  the  spectrum  used  for  this  work. 


\  Earth’s 
\  Rotation 

/ 

Atmosphere 

Figure  1.6  Diagram  of  air  mass  geometry 


The  AM  1.5  standard  spectrum  includes  all  of  the  major  effects  of  scattering  and 
absorption  that  occur  as  light  passes  through  the  Earth’s  atmosphere.  It  is  calculated  for  a 
plane  with  surface  normal  directed  at  the  Sun  and  tilted  toward  the  equator  at  37°,  chosen 
because  it  represents  an  average  latitude  for  the  48  contiguous  United  States.  It  is 
calculated  based  on  a  standard  atmosphere  for  the  United  States  with  average 
temperature,  pressure,  aerosol  density,  air  density,  and  molecular  species  density.  Also 
included  are  estimates  of  turbidity,  water  vapor  content,  ozone  content,  and  surface 
spectral  albedo  (reflectivity  of  the  soil).  The  contributions  of  some  of  these  effects  are 
labeled  in  Figure  1.5. 
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Numerical  Optimization22 

Optimization  is  the  problem  of  finding  a  global  maximum  or  minimum  of  a  given 
function.  These  two  processes  are  equivalent  because  maximizing  a  function  f(x)  is  the 
same  as  minimizing  -f(x).  There  are  many  well-known  optimization  algorithms,  but  I  use 
only  two  that  are  relatively  straightforward  for  the  present  work.  One,  called  Brent’s 
method  ,  is  useful  for  solving  one-dimensional  optimization  problems.  The  other  is 
known  as  the  simplex  method  because  it  solves  multi-dimensional  problems  by 
manipulation  of  a  geometrical  structure  called  a  simplex  along  the  topology  of  the 
function.  Both  these  methods  are  straightforward  in  that  they  do  not  require  the 
derivative  of  the  function  to  be  calculated.  Furthermore,  they  are  easily  implemented  in 
code;  in  fact,  Mathworks  includes  them  as  built-in  functions  in  Matlab. 


Figure  1 .7  Graphical  representation  of  the  golden  section  search  (left)  and  parabolic 
interpolation  (right).  A  and  D  bracket  the  search  interval,  C  is  the  current  known 
minimum.  The  next  iteration  will  try  point  B.  Since  B  is  less  than  C,  B  will 
replace  C  and  C  will  replace  D;  if  B  were  greater  than  C,  B  would  replace  A. 

Brent’s  method  uses  a  combination  of  the  golden  section  search  and  parabolic 
interpolation  to  quickly  and  robustly  find  a  minimum  with  the  fewest  possible  number  of 
function  evaluations.  The  method  requires  that  the  minimum  be  bracketed  between  two 
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points.  Golden  section  search  is  a  slow  but  sure  way  to  converge  on  the  minimum  within 
the  bracket  without  fail.  It  always  moves  from  the  current  known  minimum  into  the 
larger  segment  of  the  bracketed  interval,  as  shown  in  Figure  1.7.  The  distance  moved 

3  —  -v/5 

from  the  current  position  is  always  a  fraction  equal  to  — — — ,  which  is  known  as  the 

golden  section.  It  can  be  shown  that  using  the  golden  section  will  minimize  the  worst 
case  error  taking  the  step  might  cause. 

On  the  other  hand,  parabolic  interpolation  may  converge  on  the  minimum  in  only 
a  few  iterations.  Then  again,  it  might  flail  about  wildly  and  never  converge  at  all.  This 
procedure  involves  fitting  a  parabola  to  the  bracketing  interval  and  the  current  known 
minimum,  then  choosing  the  minimum  of  the  parabola  to  be  the  new  minimum  of  the 
function.  In  general,  there  is  no  reason  to  believe  that  the  function  can  be  well 
approximated  by  a  parabola  at  all,  which  is  why  this  method  can  fail  so  spectacularly.  Of 
course,  if  we  are  sufficiently  close  to  the  minimum  of  the  function,  or  if  the  function  is 
sufficiently  smooth,  Taylor’s  theorem  tells  us  that  a  parabola  will  be  a  good 
approximation.  The  trick  is  to  use  parabolic  interpolation  as  long  as  it  behaves,  and 
switch  to  the  golden  section  search  otherwise.  Brent’s  method  cleverly  accomplishes  this 
task  in  a  robust  and  efficient  way. 

When  the  dimensionality  of  the  problem  is  greater  than  one,  the  simplex  method 
must  be  used.  A  simplex  is  a  geometric  structure  that  exists  in  //-dimensional  space  with 
R+l  vertices.  In  two-space  it  is  a  triangle,  and  in  three-space  it  is  an  irregular 
tetrahedron.  In  optimization,  each  vertex  is  a  point  somewhere  on  the  topology  of  the 
function.  On  each  iteration,  a  new  vertex  is  chosen  for  the  simplex  that  attempts  to  move 
the  structure  downhill.  Normally,  downhill  movement  is  accomplished  by  reflection, 
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where  the  highest  vertex  is  given  up  for  a  lower  one  opposite  it,  or  contraction  toward  the 
lowest  vertex. 

Figure  1.8  shows  an  example  function  for  the  current  problem.  The  details  of 
how  it  was  generated  and  the  relevant  constraints  are  discussed  in  Chapter  2.  Suffice  to 
say  that  it  represents  a  two-dimensional  optimization  problem,  and  the  third  dimension  in 
the  figure  is  present  simply  to  emphasize  the  topology  of  the  function.  The  simplex  is  a 
triangle  in  the  two-dimensional  space,  and  each  iteration  will  (hopefully)  move  the 
triangle  closer  to  the  floor  of  the  valley.  Note  that  the  simplex  algorithm  assumes  infinite 
space,  so  I  had  to  impose  the  problem’s  constraints  artificially.  The  function  is  undefined 
outside  of  the  boundary  formed  by  the  line  with  x-  and  y-intercept  equal  to  0.5,  as  well  as 
anywhere  outside  of  the  first  quadrant.  The  steep  cliff  at  the  boundary  is  the  result.  I 
decreased  the  height  of  the  cliff  for  easier  visualization.  More  discussion  about  the 
difficulties  caused  by  artificially  imposing  the  constraints  can  be  found  in  Chapter  2. 


Figure  1.8  Graphical  representation  of  a  simplex  on  the  function  to  be  minimized. 
On  each  iteration,  a  new  vertex  will  be  chosen  for  the  simplex  that  attempts  to 
move  it  toward  the  valley  floor. 
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From  this  introduction,  it  is  easy  to  see  how  nanoshells  can  be  tuned  over  the 
optical  spectrum  for  use  in  solar  applications.  However,  it  is  not  clear  which  nanoshells 
would  be  best,  or  in  what  proportion  to  mix  them.  In  the  following  chapters  I  expand  on 
how  I  used  the  optimization  algorithms  described  above  to  find  answers  to  these 
questions. 
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Chapter  2  - 

Description  of  the  Algorithm 


The  goal  of  this  project  is  to  calculate  the  optimal  mixture  of  different  nanoshell 
species  (core  and  shell  sizes)  to  absorb  or  scatter  solar  energy.  I  accomplish  this  goal  by 
using  Mie  theory  coupled  with  standard  optimization  algorithms  as  implemented  in 
Matlab.  I  also  give  a  brief  overview  of  effective  techniques  for  using  my  code  to  solve 
such  problems  and  discuss  implicit  approximations  required  to  make  this  problem 
tractable. 


Overview  of  the  Code 

The  Matlab  code,  as  listed  in  Appendix  A,  is  liberally  commented,  so  I  will  only 
describe  it  here  briefly.  The  flow  is  as  follows: 


1)  assign  values  for  constants  and  parameters, 

2)  open  and  read  data  files  that  contain  particle  cross  sections, 

3)  calculate  the  mixture  that  maximizes  absorbed  or  scattered  light 
for  each  unique  combination  of  particle  species, 

4)  track  the  combination  of  species  that  gives  optimal  results,  and 

5)  report  the  answer. 


The  user  defines  the  desired  core  sizes,  shell  thicknesses,  shell  material  (silver  or  gold), 
and  dielectric  embedding  medium  to  be  included  in  the  optimization.  Filenames  for  Mie 
calculations  are  assembled  from  this  information  in  a  standard  format,  but  can  be  easily 
changed  to  suit  the  individual  needs  and  preferences  of  the  user. 

My  algorithm  searches  the  path  specified  by  the  variable  ‘shell_path’  for  the 
required  data  files.  If  any  data  is  missing,  the  program  creates  a  file  named 
‘required  runs.out’  with  a  list  of  missing  files  and  exits.  Otherwise,  the  program 
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maintains  the  data  files  in  memory  so  as  to  avoid  time-consuming  disk  accesses.  For  a 
well  designed  run,  the  memory  required  is  on  the  order  of  1  MB. 

All  possible  unique  combinations  of  data  are  passed  to  the  ‘findlocalmin’ 
function,  where  the  user  is  able  to  set  various  convergence  parameters  and  select  between 
Brent’s  method  for  one-dimensional  optimization  or  the  simplex  method  for  multi¬ 
dimensional  problems.  More  information  about  the  built-in  Matlab  functions  used  here 
and  the  convergence  parameters  that  control  them  can  be  found  on  the  Mathworks 
website'  .  Both  methods  are  capable  of  solving  the  one-dimensional  problem,  and  both 
arrive  at  the  same  solution,  as  expected.  However,  Brent’s  method  is  simpler  and  faster 
than  the  simplex  method. 

The  function  ‘residual’  is  minimized  by  one  of  the  above  methods.  It  calculates 
the  amount  of  light  not  absorbed  or  scattered  by  the  particles.  Written  mathematically, 
the  function  is 

R 

where  /  is  the  solar  irradiance,  c  is  the  cross  sectional  area  physically  covered  by  particle 
species  i,  q  is  the  efficiency  of  the  particle,  R  is  the  number  of  particles  in  the  mixture, 
and  the  integral  is  taken  over  wavelength.  For  each  combination  of  species,  this  equation 
is  minimized  using  the  cfs  as  the  optimization  parameters.  The  user  specifies  the  total 
coverage  of  the  particles  on  the  surface  as  a  constant  ( Zc,  =  c total  )•>  so  there  are  R-l 

R 

dimensions  or  degrees  of  freedom  in  the  problem.  The  function  checks  that  the  total 
absorption  or  scattering  cross  sections  for  all  the  particles  on  the  surface  is  less  than  unity 
( 1  -  Zc,<7,  >  0 )  to  ensure  that  a  unit  area  of  the  surface  does  not  absorb  or  scatter  more 

R 
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light  than  is  incident  upon  it.  If  the  inequality  is  false  for  a  given  wavelength,  the 
absorbed  or  scattered  light  is  capped  by  whatever  is  available  in  the  AM  1.5  spectrum. 

From  a  general  viewpoint,  the  code  is  useful  in  any  situation  that  requires  finding 
an  optimal  mixture  of  component  curves  to  approximate  a  target  curve.  For  example,  this 
work  focuses  exclusively  on  optimally  converting  solar  energy  into  other  forms,  so  the 
ASTM  AM  1.5  standard  solar  spectrum"  is  loaded  as  a  target  curve.  Any  target  curve 
whatsoever  would  be  equally  acceptable.  Additionally,  all  of  the  particle  cross  sections, 
which  are  the  component  curves,  are  calculated  externally  and  supplied  as  data  files. 

This  arrangement  allows  maximum  flexibility  in  adapting  the  code  to  future  applications. 
For  instance,  one  possibility  is  the  calculation  of  optimal  mixtures  of  asymmetrical 
nanoshells  to  fit  an  experimental  nanoshell  spectrum  as  an  alternative  to  the  much 

26  27 

debated  mean  free  path  dielectric  function  modification  ’  . 

Use  of  the  Code 

N\ 

Given  that  there  are - unique  combinations  of  particle  species,  where  N 

(N-R)!R! 

is  the  number  of  different  species,  some  finesse  is  recommended  in  selecting  species  for 
trial  so  that  N  is  minimized.  Judicious  choice  of  particle  geometry  can  dramatically  help 
the  situation.  Generally  speaking,  large  particles  tend  to  scatter  and  small  particles  tend 
to  absorb.  Therefore,  it  makes  sense  to  limit  overall  particle  size  to  the  most  logical 
regime  for  the  problem  at  hand.  A  guideline  might  be  that  absorption  problems  will 
require  particles  smaller  than  150  nm,  and  scattering  problems  require  particles  larger 
than  100  nm.  Additionally,  the  spectral  peak  position  is  partially  controlled  by  shell 
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thickness,  so  it  makes  sense  to  maximize  the  range  of  that  parameter  using  the  fewest 
core  and  shell  sizes. 

Another  technique  is  to  start  with  sparse  size  increments  and  hone  in  on  the 
correct  result  in  future  runs.  For  example,  if  I  want  to  find  the  best  single  particle  for 
absorbing  sunlight  with  a  surface  coverage  of  23%,  I  might  start  by  choosing  four  cores 
sizes  (0,  50,  100,  and  150  mn)  and  six  shell  thicknesses  (5,  10,  15,  20,  25,  and  30  mn).  If 
I  consider  both  gold  and  silver  shells,  there  are  48  unique  combinations.  Finding  that  the 
best  answer  is  [ri,  x{\  =  [50,  60]  mn,  I  might  try  core  sizes  in  5  mn  increments  from  25  mn 
to  75  mn  and  shell  thicknesses  in  1  mn  increments  from  7  to  13  mn.  After  checking  a 
total  of  only  202  unique  combinations,  I  will  find  the  best  answer  to  be 
[i'i,  n]  =  [40,  51]  mn.  This  is  almost  an  eighth  the  number  of  checks  necessary  to  search 
the  whole  space  in  5  mn  core  size  increments  and  1  mn  shell  size  increments. 

It  is  important  to  note  that  the  technique  may  settle  on  a  local  minimum, 
completely  missing  a  sharp  global  minimum  lurking  between  the  core  and  shell  sizes 
included  in  the  original  run.  If  the  spacing  between  shell  thicknesses  is  small  enough, 
that  outcome  will  usually  be  avoided;  sparsely  spaced  cores  tend  to  be  acceptable.  From 
an  engineering  standpoint,  the  computational  time  saved  must  be  weighed  against  the 
possibility  of  missing  the  optimal  solution. 

After  selecting  the  shell  species  to  run,  it  is  important  to  consider  the  total  surface 
coverage  expected  for  the  system.  Surface  coverage  is  a  major  constraint  that  can 
drastically  change  the  optimization  results.  Figure  2.1  shows  how  the  minimization 
function,  and  therefore  the  optimal  mixture  ratio  for  given  nanoshell  species,  changes  as 
total  surface  coverage  is  increased  for  a  representative  example.  There  is  no  guarantee 
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species  chosen  for  a  particular  coverage  will  continue  to  be  the  best  choice  if  coverage  is 


increased.  Without  a  specific  scenario  in  mind,  I  can  only  outline  general  guidance  about 


what  to  choose  for  this  parameter. 


[r,,  r,]  =  [25,  40]  nm,  Ag  Nanoshell 


Figure  2. 1  Variation  of  minimization  function  for  mixture  of 
silver  shells  with  [ri,  r2]  =  [25,  40]  nm  and  gold  shells 
with  [ri,  r2]  =  [55,  65]  nm  as  coverage  is  increased 


Experimentally,  we  regularly  find  interparticle  spacing  of  about  one  particle 
diameter  when  we  deposit  nanoshells  on  a  surface.  An  upper  limit  of  surface  coverage  in 
this  regime  is  easily  found  by  calculating  coverage  for  a  hexagonal  array  of  nanoparticles 
with  spacing  of  one  diameter.  The  result  is  coverage  of  approximately  23%,  as  shown  in 
Appendix  B.  We  also  have  the  ability  to  create  close-packed  arrays  of  particles  that 
cover  up  to  62%  of  the  surface  .  Note  that  high  surface  coverage  can  reduce  the 
accuracy  of  the  result  because  of  particle  coupling  issues  as  discussed  below. 

Finally,  consider  the  choice  of  convergence  parameters  used  in  the 
‘findlocalmin’  function.  Brent’s  algorithm  tends  to  converge  regularly  regardless  of 
the  parameters  chosen,  but  the  simplex  method  can  be  more  finicky.  Part  of  the  problem 
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is  the  sudden  change  in  the  residual  at  the  boundary  of  the  region  where  the  function  to 
be  minimized  is  defined.  Brent’s  method  poses  no  problem  because  the  minimum  is 
bracketed  ahead  of  time,  but  the  simplex  method  only  recognizes  an  infinitely  continuous 
space.  That  means  that  it  will  try  solutions  that  are  not  allowed  by  the  surface  coverage 
constraints.  For  example,  if  the  total  surface  coverage  is  specified  to  be  ctota/  =  0.3,  it 
may  try  a  solution  with  a  =  0.4  and  c?  =  -0.1.  The  code  prevents  unphysical  results  by 
automatically  setting  the  residual  to  the  maximum  value  (no  light  is  absorbed  or 
scattered)  when  the  optimization  values  fall  outside  the  permissible  domain. 

When  the  minimum  is  located  on  this  boundary,  the  simplex  method  does  not 
converge  quickly.  It  is  possible  to  solve  this  problem  by  increasing  the  maximum 
allowable  number  of  iterations  and  function  evaluations,  but  unbearably  long  execution 
times  will  result.  A  better  solution  is  to  increase  the  termination  tolerance,  ‘TolX,’  which 
controls  the  accuracy  the  computer  requires  for  a  solution.  Usually,  convergence 
continues  until  the  answer  is  within  the  square  root  of  the  floating  point  accuracy  of  the 
machine,  as  that  value  is  the  best  the  computer  can  accomplish"  .  In  this  case,  we  are 
optimizing  for  mixtures  of  nanoshells  with  poorly  known  concentration  in  solution,  so  it 
makes  no  sense  to  find  optimal  surface  coverage  to  eight  decimal  places.  That  precision 
would  be  impossible  to  reproduce  when  making  real  samples  in  the  laboratory. 

Inherent  Approximations 

As  stated  previously,  particle  cross  sections  for  each  shell  species  are  calculated 
externally  and  supplied  as  data  files.  I  employed  Mie  scattering  theory  to  calculate  the 
particle’s  absorption  and  scattering  efficiencies  as  a  function  of  wavelength.  This  theory 
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does  an  excellent  job  of  predicting  the  scattering  and  absorption  cross  sections  of  dilute 
nanoshell  suspensions  in  a  homogenous  mediunr  .  Unfortunately,  there  are  several 
confounding  factors  that  make  it  difficult  to  naively  apply  Mie  theory  to  the  current 
optimization  problem. 

First  is  the  fact  that  many  energy  conversion  problems,  particularly  the  conversion 
of  electromagnetic  to  electrical  energy,  require  the  plasmonic  particles  to  be  deposited  on 
a  surface.  Naturally,  the  homogeneity  of  the  environment  is  destroyed  by  the  presence  of 
the  surface.  The  asymmetry  of  the  surrounding  medium  causes  an  increase  in  the  net 
dipole  moment  of  higher  order  plasmon  resonance  modes,  which  allows  light  to  couple 
with  those  modes  more  efficiently.  This  effect  is  particularly  strong  for  the  component  of 
the  incident  electric  field  oriented  perpendicular  to  the  surface'  .  As  a  result,  Mie  theory 
may  underestimate  the  strength  of  the  quadrupole  mode  relative  to  the  dipole  in  a 
complicated  way  that  depends  on  the  incidence  angle  and  polarization  of  the  light. 
Furthermore,  the  resonance  energy  will  shift  because  of  hybridization  with  surface 
plasmon  states  if  the  surface  has  free  electrons  in  its  conduction  band  . 

A  crude,  but  useful,  approximation  of  these  effects  is  to  assume  a  homogenous 
effective  medium6  for  the  particle  environment  by  calculating  a  weighted  average  of  the 
surface  and  medium  dielectric  functions  ( £eff(co )  =  wl£surf(co)  +  w2smed).  Appropriate 

weights  can  be  determined  geometrically  if  the  particle’s  optical  properties  are  affected 
only  by  the  dielectric  environment  within  one  particle  radius  .  The  weights  would  then 
be  equal  to  the  ratio  of  the  volume  of  surface  to  the  total  volume  within  that  distance  and 
the  volume  of  medium  to  the  total.  Figure  2.2  shows  the  geometry  involved.  Although 
this  solution  is  less  than  ideal  because  it  cannot  accurately  model  the  affect  of  the  surface 
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on  relative  peak  height  and  shape,  it  does  partially  correct  for  changes  in  peak  position. 
This  design  tradeoff  enables  quick  calculation  of  absorption  and  scattering  efficiencies 
for  extremely  large  numbers  of  particles  with  standard  Mie  theory,  at  the  expense  of 
modeling  the  quadrupole  peak  well. 


The  volume  of  the  surface  that  is  within 
1  radius  of  the  particle  is 

V  =  j  ml  (6r2  -  ^2 )  =  f  m'i 
Therefore, 


1  -  r2 )  28 

and 


w2  =  1  -  Wj 


23 
28  ' 


Figure  2.2  Dielectric  environment  of  a  nanoshell  on  a  surface  -  the  environment 
within  1  particle  radius  contributes  to  the  optical  response 


The  second  problem  is  the  close  proximity  of  particles  on  the  surface.  The 
particles  cannot  necessarily  be  treated  as  independent.  It  is  clear  from  conservation  of 
energy  that  when  two  particles  are  close  enough  for  their  extinction  cross  sections  to 
overlap,  light  incident  on  the  overlapping  area  is  not  extinguished  by  both  particles.  That 
is,  it  cannot  be  absorbed  twice.  Therefore,  the  total  extinction  cross  section  of  the  two  is 
not  simply  the  sum  of  the  independent  cross  sections.  Hybridization  between  the 
particles  via  columbic  forces  can  change  the  scattering  profile  so  completely  that  it  is 
insufficient  to  correct  the  cross  section  by  subtracting  the  overlap  area  from  the  total34. 

oc  -7/: 

Furthermore,  conducting  substrates  can  enhance  coupling  between  particles  ’  . 

In  order  to  gain  some  insight  about  the  magnitude  of  the  error  introduced  by 
plasmonic  coupling,  I  investigated  these  effects  for  a  specific  example  using  the  T-matrix 
method  with  code  that  is  freely  available  from  NASA  .  Figure  2.3  shows  the  probability 
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light  will  scatter  from  a  dimer  system  as  a  function  of  interparticle  spacing.  The  distance 
dependence  is  oscillatory  so  that  sometimes  the  probability  is  enhanced  and  sometimes  it 
is  diminished.  In  a  worst  case  scenario,  assuming  independent  particles  will  result  in  a 
misestimate  of  the  scattering  probability  by  4%. 


Figure  2.3  Comparison  of  the  probability  a  gold  colloid  dimer  will  scatter  light  with 
558  nm  wavelength  if  the  particles  are  treated  as  independent  scatterers  versus  if 
they  are  treated  as  a  coupled  two  particle  system;  the  particles  have  radii  equal  to 
82  nm  and  60  nm 

To  estimate  the  number  of  particles  that  will  be  affected  by  coupling,  consider 
that  interparticle  spacing  will  be  governed  by  Poisson  statistics,  as  shown  in  Appendix  C. 
Therefore,  the  most  probable  number  of  particle  couplings  in  a  unit  area  is 

1  2  j  2 

a  , 

where  x  is  the  concentration  of  particles  on  the  surface  and  d  is  a  parameter  that  defines 
the  interparticle  spacing  at  the  limit  of  the  validity  of  independent  scattering.  The 
difficulty  of  this  problem  lies  in  the  determination  of  the  parameter  d  and  the  efficient 
calculation  of  dimer  cross  sections  for  the  statistical  distribution  of  interparticle  distances 
less  than  d.  Spacing  greater  than  two  particle  diameters  is  a  rough  estimate  of  the 
distance  required  for  interparticle  coupling  effects  to  be  negligible  .  The  best  way  to 
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reduce  the  error  introduced  is  to  keep  the  nanoparticle  surface  coverage  low,  although 
results  are  still  meaningful  for  higher  coverage  if  one  is  mindful  of  the  approximations 
involved. 

The  last  complicating  factor  is  multiple  scattering  effects.  This  factor  refers  to  the 
possibility  that  light  scattered  by  one  particle  is  subsequently  absorbed  by  another. 
Fortunately,  consideration  of  the  angular  dependence  of  the  differential  scattering  cross 
section  reveals  that  scattering  occurs  predominately  in  the  forward  and  backward 
directions,  as  shown  in  Figure  2.4.  Since  the  particles  that  might  subsequently  absorb 
scattered  light  lie  at  right  angles  to  the  direction  of  incidence,  I  expect  multiple  scattering 
to  be  a  negligible  effect.  Furthermore,  the  effect  will  be  an  issue  only  for  large 
nanoshells  because  smaller  particles  tend  to  scatter  less. 


Dipole 

Quadrupole 


Figure  2.4  Angular  dependence  of  scattering  for  unpolarized  light  incident 
normal  to  the  surface  with  wavelength  equal  to  the  dipole  resonance  and 
the  quadrupole  resonance;  a  homogenous  effective  medium  is  assumed 


The  algorithm  described  here  uses  various  approximations  that  are  necessary  to 
make  computation  time  realistic.  While  accuracy  suffers  slightly  as  a  result,  the  solution 
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still  contains  considerable  information.  The  quadrupole  mode  in  particular  will  be 
susceptible  to  error  because  of  these  approximations.  However,  it  is  unlikely  that  a 
ideally  optimal  solution  will  differ  significantly  from  the  reported  results. 
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Chapter  3  - 

Two  Scenarios  in  Detail 

I  used  my  optimization  algorithm  to  solve  two  scenarios  of  practical  interest.  The 
first  scenario  is  using  plasmonic  particles  to  couple  sunlight  into  a  silicon  photodiode, 
and  the  second  is  to  maximize  visible  light  scattering  from  a  glass  surface.  These  are  two 
of  many  interesting  applications  and  are  only  meant  to  serve  as  an  example  of  the 
possibilities. 

Absorption  of  Sunlight  on  a  Silicon  Photodiode 

In  2005,  Schaadt  et  al.  demonstrated  that  gold  colloid  deposited  on  a  silicon  p/n 
junction  can  enhance  the  photocurrent  by  80%  with  surface  coverage  as  low  as  1%.  The 
mechanism  for  this  effect  is  not  well  understood,  but  it  is  clear  from  unpublished 
experiments  that  it  only  occurs  for  absorptive  particles.  If  colloid  can  couple  light  into 
the  junction  so  dramatically,  it  should  be  possible  to  take  advantage  of  the  tunability  of 
nanoshells  to  maximize  the  absorption  of  the  incident  light  and  fully  exploit  the 
enhancement. 

My  algorithm  searched  for  the  best  mixture  of  nanoshells  for  absorbing  AM  1 .5 
sunlight  with  wavelengths  between  300  and  1100  mn.  Wavelengths  longer  than  1100  mn 
have  energy  less  than  the  band  gap  of  silicon  and  will  not  participate  in  the  generation  of 
photocurrent.  Both  gold  and  silver  particles  were  allowed.  I  used  Johnson  and 

TO 

Christy’s  data  for  the  dielectric  function  of  gold,  the  CRC  Handbook  of  Chemistry  and 
Physics40  for  silver,  and  the  Handbook  of  Optical  Constants  of  Solids41  for  silicon.  Since 
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small  plasmonic  particles  tend  to  be  the  best  absorbers,  I  focused  on  ri  <  85  nm,  and 
surface  coverage  was  capped  at  40%. 

The  optimal  result  is  a  mixture  of  only  two  species,  both  gold  nanoshells,  that 
absorbs  582  W/m2  out  of  805  W/m2  available.  The  first  species  is  [ri,  r2]  =  [47,  58]  nm 
with  coverage  of  27.8%,  and  the  other  is  [ri,  r2]  =  [28,  42]  nm  with  coverage  of  12.2%. 
Correcting  for  particle  size,  it  will  require  six  of  the  larger  species  for  every  five  of  the 
smaller  to  achieve  the  desired  coverage  ratio.  Mixing  the  species  in  this  6:5  ratio  will 
enable  absorption  of  nearly  three-quarters  of  the  incident  irradiation  as  shown  in  Figure 
3.1.  Allowing  mixtures  of  three  or  more  species  does  not  increase  the  absorbed  energy 
unless  surface  coverage  is  also  allowed  to  increase.  If  total  surface  coverage  is  increased 
to  50%,  an  additional  96  W/m2  can  be  harvested  from  the  light  as  shown  in  Figure  3.2. 

On  a  real  surface,  interparticle  effects  would  enable  the  particles  to  absorb  the  low  energy 
light  that  is  considered  “residual”  by  my  algorithm. 


-AM  1.5  Solar  Spectrum 
-  Absorption  of  NS  Mixture 
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Figure  3.1  Optimal  spectra  of  light  absorbed  for  a  silicon  surface  40%  covered  by 
particles 
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-AM  1.5  Solar  Spectrum 
Total  Absorption  of  Mixture 
'  [r,,  rj  =  [47,  58]  nm,  Au  32.2% 


-  [rl,  r2]  =  [47,  58]  nm,  Au  32.2% 

-  [rl,  r2]  =  [28,  42]  nm,  Au  9.3% 
-r  =  32,  Au  Colloid  8.5% 


Figure  3.2  Optimal  spectra  of  light  absorbed  for  a  silicon  surface  50%  covered  by 
particles 

It  is  interesting  to  consider  how  these  “black”  surfaces  compare  to  a  surface 
painted  black  with  some  traditional  surface  coating.  Obviously,  black  paint  is  nearly  a 
perfect  absorber  of  light  in  the  visible  wavelength  region.  In  order  to  keep  the 
comparison  fair,  I  will  consider  a  submonolayer  of  paint  particles  with  total  surface 
coverage  equal  to  40%.  Black  pigment  is  almost  always  made  from  carbon  ,  and  Mie 
scattering  theory  was  successfully  used  to  describe  clusters  of  soot  by  many 
researchers4’.  I  will  model  black  paint  as  small  (20  nm)  spherical  particles  and  use  the 
complex  index  of  refraction  of  soot  (1.75+i0.435)  reported  by  d’Almeida  et  al.44.  Surely, 
this  leap  can  be  no  worse  than  modeling  the  radar  cross  section  of  sparrows  using  bird¬ 
sized  balls  of  water45. 

Figure  3.3  shows  the  absorption  efficiency  of  carbon  soot  from  300  to  1100  nm. 
The  optical  behavior  is  relatively  flat  across  the  wavelength  region  of  interest  with  a 
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distinct  lack  of  resonance  behavior.  The  black  color  of  soot  originates  from  this  even  and 
smooth  absorption  spectrum.  However,  the  optical  efficiency  is  extremely  low.  As  a 
result,  at  40%  surface  coverage,  the  particles  only  absorb  a  miniscule  48  W  out  of  805  W 
available  in  a  square  meter  of  sunlight.  Some  studies  show  that  absorption  cross  sections 
of  soot  aggregates  are  vastly  enhanced  by  interparticle  effects43,  and  in  absolute  terms, 
thick  coats  of  black  paint  have  absorptivity  greater  than  90%46.  But,  in  terms  of 
absorbing  solar  energy  with  sparse  surface  coverage,  black  paint  simply  cannot  compete 
with  plasmonic  particles. 


- AM  1.5  Solar  Spectrum 


Figure  3.3  Absorption  spectrum  of  soot 


Scattering  of  Sunlight  from  a  Glass  Surface 

The  second  scenario  I  will  discuss  is  optimizing  light  scattering  from  a  glass 
surface.  This  problem  is  interesting  in  that  it  is  the  inverse  of  the  absorption  problem, 
and  because  it  may  find  useful  application  in  window  treatments  or  other  surface 
modifiers.  The  Romans  used  the  attractive  optical  scattering  properties  of  gold  and  silver 
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colloid  in  the  4th  century  for  coloring  glass.  One  of  the  most  famous  examples  is  the 
Lycurgus  Cup,  on  display  at  the  British  Museum.  The  cup  normally  appears  green,  but 
changes  color  to  red  if  viewed  with  a  light  source  in  transmission. 

In  this  example,  I  optimize  for  wavelengths  limited  to  the  visible,  between  380 
and  820  nm.  I  took  the  dielectric  constant  of  glass  to  be  2.04  over  the  entire  domain  of 
wavelengths.  As  before,  both  gold  and  silver  particles  are  allowed,  and  their  complex 
index  of  refraction  was  taken  from  the  same  sources  used  previously.  I  expected  large 
particles  to  be  the  optimal  scatterers,  so  I  concentrated  on  nanoshells  with  cores  larger 
than  160  nm  diameter,  and  large  colloid. 


-AM  1.5  Solar  Spectrum 
-  Scattered  Light 


-  Scattering 


Figure  3.4  Optimal  spectra  of  light  scattered  for  a  glass  surface  30%  covered  by 
particles 


Initially,  I  permitted  total  surface  coverage  to  be  40%,  but  it  turns  out  that 

2 

coverage  of  30%  is  sufficient  to  scatter  almost  all  of  the  incident  irradiation:  571.5  W/m 
out  of  578  W/m  .  I  was  further  surprised  to  find  that  a  single  particle  species  is  all  that  is 
required  to  achieve  this  feat.  The  best  particle  for  the  job  is  silver  colloid  with  a  diameter 
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of  210  nm.  Its  scattering  efficiency  and  spectrum  under  AM  1.5  irradiation  is  shown  in 
Figure  3.4. 


- AM  1.5  Solar  Spectrum 

- Scattered  Light 
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Figure  3.5  Scattering  spectrum  of  TiCF 

To  gain  some  insight  about  how  the  performance  of  this  particle  compares  to  the 
visible  light  scattering  of  everyday  objects,  I  show  the  scattering  spectrum  of  the  main 
ingredient  of  white  paint.  White  pigment  can  be  made  from  a  variety  of  oxides,  but  one 
of  the  best  is  titanium  dioxide42.  Unlike  carbon,  titanium  dioxide  has  very  high  extinction 
efficiency  with  pronounced  peaks,  as  shown  in  Figure  3.5.  I  used  the  index  of  refraction 
listed  in  the  Flandbook  of  Optical  Constants  of  Solids47.  Titanium  dioxide  does  not 
scatter  the  low  energy  light  well,  and  with  30%  surface  coverage  it  only  scatters 
458  W/m".  However,  TiC>2  will  still  appear  to  be  more  brilliantly  white  than  silver 
colloid  because  its  absorption  spectrum  is  nearly  zero  for  all  visible  wavelengths.  Silver 
colloid  absorbs  blue  light,  making  it  look  reddish-brown.  The  fact  that  TiC>2  does  not 
scatter  red  light  well  may  leave  it  with  a  blue  tinge,  but  silver  will  still  be  less  white.  The 


-  Scattering 

-  Absorption 
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code  does  not  correct  for  the  absorption  spectrum  when  optimizing  for  scattering. 
Including  interparticle  effects  would  eliminate  the  difficulties  caused  by  that  omission. 

In  summary,  both  example  scenarios  show  that  plasmonic  nanoparticles  are 
excellent  tools  for  designing  optical  behavior.  They  can  be  used  in  both  absorption  and 
scattering  applications  for  converting  solar  energy  into  other  usable  forms.  Gold 
nanoshell  species  [ri?  r2]  =  [47,  58]  nm  mixed  with  [ri,  r2]  =  [28,  42]  nm  is  optimal  for 
absorbing  sunlight,  and  silver  colloid  with  r  =  105  nm  is  best  for  scattering.  The 
particles’  naturally  broad  ensemble  resonances  are  ideal  for  interacting  with  broadband 
sunlight.  Using  this  optimization  code  allows  researchers  to  quickly  detennine  the  best 
particle  for  a  given  experiment,  and  hopefully  it  will  enable  informed  design  choices  for 
nanotechnological  optical  devices  of  the  future. 
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Chapter  4  - 

Conclusions  and  Future  Work 

Physical  understanding  of  how  nanoshells  interact  with  their  environment  could 
be  used  to  improve  the  algorithm  as  it  now  stands.  Some  examples  of  the  physical 
processes  involved  are  discussed  in  Chapter  2.  Future  improvements  to  the  code  should 
incorporate  those  effects.  Additionally,  an  engineering  improvement  that  could  decrease 
computation  dramatically  would  be  to  use  an  optimization  algorithm  to  search  all  realistic 
nanoshell  species.  It  is  subtle,  but  there  are  actually  two  optimization  problems  being 
solved.  The  first  is  optimizing  the  mixture  ratio  given  a  set  of  nanoshell  species.  The 
second  is  to  optimize  the  set  of  species.  In  the  current  algorithm,  brute  force  is  used  to 
determine  the  optimal  group  of  nanoshell  species.  If  an  optimization  algorithm  were 
substituted  for  this  brute  force  method,  the  computational  time  saved  could  be  applied  to 
more  comprehensive  calculations  of  interparticle  effects. 

In  this  work,  I  demonstrated  how  well-known  optimization  algorithms  can  be 
applied  to  solve  an  engineering  problem  in  one  nanoscale  system  -  harnessing  sunlight 
using  mixtures  of  nanoshells  on  a  surface.  As  we  learn  more  about  how  to  predict  and 
control  the  unique  properties  of  nanoshells,  applying  engineering  techniques  to  realize 
novel  devices  becomes  important.  The  code  written  for  this  project  should  be  useful  for  a 
wide  range  of  problems  as  we  attempt  to  find  new  ways  to  use  nanoshells  in  devices. 
Developing  nanoshell  applications,  and  more  generally  advancing  the  field  of 
nanotechnology,  are  truly  multi-disciplinary  endeavors:  chemistry  controls  the  properties 
of  nanoscale  matter,  physics  is  used  to  understand  them,  and  engineering  exploits  them. 
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Appendix  A 
Code 


solarfit.m 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


4; 

1700  1702 
10:55]  ; 

25  75]  ; 


num  species  = 
lambda  =  [300 
cores  =  [0  25 
shells  =  [5:5 
metal_types  = 
terms  =  3; 
coverage  =  0.4; 
interaction  =  'absb' 
medium  =  'Glass'; 


%Number  of  different  shell  species  in  mix 
1705:5:1800].';  %Wavelens,  col  vector  only 
%core  radii  to  try 
%shell  thicknesses  to  try 
['Ag';'Au'];  %shell  metal  types  to  try 
%maximum  1  quantum  number  to  use 
^fractional  unit  area  covered  by  shells 
luse  'scat'  or  'absb' 

^embedding  medium  tag 


[m_rows,m_cols]  =  size (metal_types) ; 
num  cores  =  length (cores) ;  %get  number  of 
num  shells  =  length ( shells ) ;  %get  number  of 
species  =  m  rows . *num  cores. *num  shells;  %# 


cores 
shells 
of  species 


to  try 


combinations  =  prod ([ species-  ...  %N ! / ( (N-R) ! R! )  possible  combos 
num_species+l : species] ) . /factorial (num^species) ; 

runs_needed  =  0; 

shell  path  =  ...  %path  to  Mie  calculation  data 

' /home/ j  oecole/ solar_spectrum_f itting/ aashell_data/broad/ ' ; 


load_spec_aml_5 ; 
I_total  =  sum(I_tilt) 

open  shell  files; 


%load  the  solar  spectrum 

%integrate  solar  spec  over  wavelength 

%load  all  Mie  data  into  memory 


fprintf ( 1 ,' Checking  %d  combinations  of  %d:\n' 
combinations , num  species); 


ind2  =  [ ] ; 
count  =  0; 

resid  light  min  =  I  total; 
for  ind  =  [ 1 : combinations ] 
if (mod (ind, 200) ==0) 
fprintf  (1,  '  .  '  ) 
count  =  count  +  1; 
if (count==70) ,  fprintf (1, 

end 


^initialize  to  max  possible  value 

%try  every  combination 

%give  user  screen  output 

iso  they  know  something' s  happening 


\n '  )  ,  count=0;,  end 


ind2  =  next  combo ( species , num  species , ind2 ) ; 

%get  indices  for  Mie  data  to  use  in  this  iteration 
[ad2, resid_light (ind) ]  =  ... 

find  local  min (I  tilt, coverage, C  ef f ( : , ind2 ) ) ; 
%find  optimal  coverage  ratio  for  this  mixture 

if (resid  light (ind)  <=  resid  light  min) 

%if  this  is  the  best  mix  so  far  . . . 

resid  light  min  =  resid  light (ind);  %save  the  answer 
ad2  min  =  ad2 ;  %save  the  optimization  parameters 

ind  min  =  ind2;  %save  the  indices 

end 
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next  iteration. . . 


52  end  % 

53 

54  fprintf ( 1 , ' Done . \n ' ) 

55 

56  fprintf (1, ' \nBest  answer  is:\n')  %write  the  best  answer  to  screen 

57  for  m=[l:num  species] 

58  [i,j,k]  =  ind2 sub ( [myrows, num  cores, num  shells], ind  min (m) ) ; 

59  fprintf (1, ' \tRl  =  %3d,  R2  =  %3d.  Shell  Type  =  %s.  Coverage  = 

60  %5.2f%%\n', 

61  cores (j) ,shells(k)+cores(j) , metal_types ( i , : ) , coverage .*100) ; 

62  end 
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find  local  min.m 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


function  [d,y]  =  find  local  min ( I , coverage, eff) 

%function  [d,y,err]  =  find  local  min ( I , coverage, eff ) 

%finds  the  optimal  density  for  mixtures  of  shell  species  on  a 
%surf ace 

%I  -  column  vector  with  elements  representing  the  irradiance  at 
%  each  wavelength 

%coverage  -  scalar  representing  fraction  of  unit  area  physically 
%covered  by  shells 
%  (pi/ (8*sqrt (3) )  is  good) 

%eff  -  matrix  of  effective  cross  sections  (num  of  wavelengths  rows 
%and  number 

%  of  shell  species  columns 

[wavelen, numspec]  =  size (eff);  %get  #  of  wavelengths  and  R 

frac  cov  =  ones ( 1 , numspec-1 ). ‘coverage . /numspec; 

%initialize  with  equal  coverage  for  all  species 

options  =  optimset ( ' Display ', 'Notify', ' TolX ' ,  . . . 

le-4  ,  ' MaxFunEvals ' ,  6000 ,  ' Max I ter ' ,  4000 ) ; 

%Set  the  convergence  parameters 

%[d,y,err]  =  fminbnd(@(x)  sum ( residual (x, coverage, I , eff )) ,  ... 

%  0, coverage, options) ; 

%  Brent's  method;  uncomment  for  1-d 

[d,y,err]  =  fminsearch ( @ (x)  sum (residual (x, coverage, I , eff )) ,  ... 

frac_cov, options)  ; 

%  Simplex  method;  uncomment  for  multi-d 
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residual.m 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


function  y  =  residual (frac_coverage, coverage, I, eff) 

%function  y  =  residual (frac_coverage, coverage, I, eff ) 

%finds  the  optimal  density  ratio  for  species  of  shells  on  a  surface 
%frac  coverage  -  Each  column  is  the  fraction  of  a  unit  area 
physically  covered  by 

%  that  shell  species.  This  is  a  row  vector  with  one  fewer 

columns 

%  than  the  number  of  columns  in  eff  or  extinct. 

%I  -  column  vector  with  elements  representing  the  irradiance  at 
%  each  wavelength 

%coverage  -  scalar  representing  fraction  of  unit  area  covered  by 
shells 

%eff  -  matrix  of  effective  scattering  or  absorbtion  cross  sections 
%  (num  of  wavelengths  rows  and  number  of  shell  species  columns) 
%The  output  is  a  column  vector  of  the  residual  light  for  each 
wavelength 

resid  =  1; 

[wavelens , num]  =  size (eff); 

frac  coverage (num)  =  coverage  -  sum(frac  coverage); 

%calculate  coverage  for  last  species  using  constraint  equation 

if ( ( f rac_coverage>=0 )  &  (frac_coverage<=coverage) ) 

%if  the  optimization  parameter  for  each  species  is  in  the  allowed 
%domain 

frac  cov  =  ones (wavelens, 1) *frac  coverage;  %expand  to  allow 

%multiply 

resid  =  1-sum ( f rac_cov . *ef f , 2 ) ;  %multiply  and  sum  on  per 

%wavelength  basis 

ind  =  find (resid  <  0);  %look  for  unphysical  results 
resid (ind)  =  0;  %fix  unphysical  results 

end 

y  =  I.*resid;  %multiply  by  solar  spectrum 
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next  combo.m 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


function  [combo]  =  next  combo (N,  R,  current  combo) 

%function  [combo]  =  next  combo (N,  R,  current  combo) 

%Calculates  the  indices  of  the  next  set  of  species 
%N  -  number  of  different  species  (total) 

%R  -  number  of  species  in  a  set 

%current  combo  -  indices  of  current  set  of  species 

if (~isempty (current  combo))  %if  already  initialized... 

if (R  ~=  length (current  combo))  %if  programming  error... 

disp (' Argument  current  combo  too  long  in  next^combo ( ) ') 
disp('<CTRL>  -  C  to  quit'); 
combo  =  zeros (1,R); 
pause 

end 

for  i  =  [ 1 : R] 

if ( (current  combo (i)  <  i)  |  (current  combo (i)  >  N-R+i)) 
%if  programming  error... 

disp (' Argument  current  combo  out-of-range  in  ... 

next^combo ( ) '  )  ; 

disp('<CTRL>  -  C  to  quit'); 

combo  =  zeros (1,R); 

pause 

end 

end 

combo  =  current  combo; 

combo (R)  =  combo (R)  +  1;  %counts  to  next  index 

for  i  =  [R:— 1:2]  %if  overflow  then  carry  the  one... 
if (combo (i)  >  (N-R+i)) 

combo ( i — 1 )  =  combo ( i — 1 )  +  1; 
for  j  =  [i:R] 

combo (j)  =  combo (j-1)  +  1; 

end 

end 

end 

if (combo (1)  >  (N-R+l))  %if  programming  error  ... 

disp (' Function  next^combo ( )  called  too  many  times.'); 
disp('<CTRL>  -  C  to  quit'); 
combo  =  zeros (1,R); 
pause 

end 

else  %if  not  already  initialized... 

combo  =  [1:R];  %initialize  for  1st  time  through  loop 

end 
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load  spec  ami  5.m 

1  load  ' ASTMG173 . csv '  -ascii  %load  standard  am  1.5  spectra 

2  lambda  indexes  =  pick  lambdas (ASTMG173 (:, 1) , lambda) ; 

3  %choose  wavelengths 

4  Etr  =  ASTMG173 (lambda  indexes, 2); 

5  %Extraterrestial  irradiance  (AMO) 

6  I  tilt  =  ASTMG173 (lambda  indexes, 3); 

7  %Irradiance  for  tilted  plane 

8  %I  direct  =  ASTMG173  2 (lambda  indexes, 4); 

9  %Direct  irradiance 
10 

11  clear  ASTMG173  lambda  indexes;  %clear  unneeded  variables 
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load  shell  data.m 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


function  [xsec]  =  load  shell  data (filename, lambdas, xsec  type) 
%function  xsec  =  load  shell  data (filename, lambdas, xsec  type) 

%  opens  a  shell  *.ext  file  and  returns  the  effective  absorption  or 
%  scattering  cross  section  for  the  shell 
%filename  -  full  path  to  needed  file 
%lambdas  -  wavelengths  needed  for  run 

%xsec_type  -  type  of  interaction  (absorption  or  scattering)  to  solve 
%  for 


if (nargin==2 )  %figure  out  which  column  is  needed 

xsec_type  =  'absb'; 

end 

if (xsec_type  ==  'scat') 
ind  =2; 

elseif (xsec_type  ==  'absb') 
ind  =  3; 

elseif (xsec_type  ==  ' extn ' ) 
ind  =  4; 

else 

disp('I  do  not  understand  the  third  argument  you  provided  to 

load_shell_data ( ) . ' ) ; 

disp('<CTRL>  -  C  to  quit'); 
pause 

end 

shell_data  =  load (filename, ' -ascii ') ; 

%load  the  required  data 

lambda  indexes  =  pick  lambdas ( shell  data, lambdas ) ; 

%select  required  wavelengths 

xsec  =  shell  data (lambda  indexes , ind) ; 

%select  required  column 
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pick  lambdas.m 

function  inds  =  pick  lambdas (data, lambda) 

%function  inds  =  pick  lambdas (data, lambda) 

%checks  data  to  make  sure  the  wavelengths  needed  for  the  run  exist 
%and  selects  the  indices  of  those  wavelengths 

%these  must  be  in  increasing  order,  because  the  function  only 
%searches  the  portion  of  data  that  hasn't  been  checked  before  in 
%order  to  shorten  the  run  time 

O, 

o 

%data  -  vector  of  wavelengths  present  in  the  data 
%lambda  -  vector  of  wavelengths  needed  for  run 

i  =  1; 

e  =  length (data) ; 
inds  =  [ ] ; 


for  1 


pick_ 


=  lambda . 1 

ind  =  find(l  ==  data(i:e))  +  i  -  1; 

%  find  next  required  wavelen 
if ( length ( ind)  ~=  1) 

disp(' Unable  to  find  specified  lambda  in  function 
lambdas ( ) ' ) ; 

disp ( ' <CTRL>-C  to  quit'); 
pause 

end 

i  =  ind+1; 

inds  =  [inds;ind];  %append  new  index  to  end  of  list 


end 
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Appendix  B- 

Surface  Coverage  Calculation 


Figure  B.  1  Hexagonal  array  of  nanoparticles 
with  one  diameter  spacing  between  them. 


The  surface  can  be  divided  into  30-60-90  triangles  with  a  base 


b  =  2r 


2  5 


and  height 


h=2^3r7. 


The  area  of  the  circle  inside  the  triangle  is 
therefore,  the  surface  coverage  is: 


—  7T 

4  n  '2  ’ 


—  7T  F  2 
4  '2 

\bh 


7T  r2 

8  a/3 /' 
n 
8a/3 


0.23. 
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Appendix  C- 

Number  of  Interparticle  Couplings  on  a  Surface 


Imagine  a  surface  with  area  A.  Particles  are  added  at  random  to  the  surface  one-by-one. 
When  the  Nth  particle  is  added,  it  causes  N-l  Bernoulli  trials.  We  define  success  to  mean 
that  the  particle  lands  within  a  distance  d  of  another  particle,  close  enough  for 
interparticle  coupling  effects  to  be  relevant.  The  probability  of  success  is 

n  d2 
P=- 


A 


The  probability  of  failure  is 


4  =  1  — 


n  d7 
~A 


As  N  increases,  the  number  of  trials  per  unit  area  is 


n  1A. 

—  =  —  >  /  = 

a  Atr 


N2  +  N 
2  A 


N 
2 A 


for  large  N.  Since  n/A  tends  toward  infinity  as  0(N~ ),  Poisson’s  approximation  will  give 
good  estimates  for  this  binomial  distribution.  The  mean  of  Poisson’s  distribution  (per 
unit  area)  is  known  to  be 

N2nd2 
A  2  A2  ’ 

which  also  equals  its  standard  distribution.  The  concentration  of  particles  per  unit  area,  x, 
can  be  written  as 

N 

leaving  us  with  the  most  probable  number  of  interparticle  couplings  per  unit  area: 

/t  =  \nx2d 2 . 

Note  that  this  is  not  the  number  of  coupled  particles  because  each  particle  can  be  coupled 
with  multiple  others.  This  estimate  includes  all  the  couplings.  We  can  find  a  lower 
bound  on  the  number  of  coupled  particles  by  noting  that  there  are 

k\ 

(  k- 2)121 

possible  two-particle  couplings  for  k  particles.  Solving  for  k,  we  arrive  at  the  quadratic 
equation: 

W^-  =  \xx2d2 
2  2 

k2  -  k-nx2d 2  =  0 

1  ±  Jl  +  47rx2d2 

- - - =  k 
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